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interacting atoms 



O 

o 

(N 

a 



Vladimir Gritsev 1 , Anatoli Polkovnikov 2 , Eugene Demler 1 
1 Department of Physics, Harvard University, Cambridge, MA 02138 
2 Department of Physics, Boston University, Boston, MA 02215 
(Dated: February 4, 2008) 

We use quantum sine-Gordon model to describe the low energy dynamics of a pair of coupled one- 
dimensional condensates of interacting atoms. We show that the nontrivial excitation spectrum of 
the quantum sine-Gordon model, which includes soliton and breather excitations, can be observed 
in experiments with time-dependent modulation of the tunneling amplitude, potential difference 
between condensates, or phase of tunneling amplitude. We use the form-factor approach to compute 
structure factors corresponding to all three types of perturbations. 
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I. INTRODUCTION 

When discussing two dimensional classical or one di- 
mensional quantum systems, no model has a wider range 
of applications than the sine-Gordon model. The sine- 
Gordon model (SG) was originally studied in the context 
of high-energy physics^ and later became a prototypi- 
cal model of low-dimensional condensed matter systems. 
In statistical physics it has been successfully applied to 
describe the Kosterlitz-Thouless transition in two dimen- 
sional superconductors and superfluids 3 -^ as well as com- 
mensurate to incommensurate transitions^ on surfaces. 
The quantum sine-Gordon model was used to describe 
the superfluid to Mott transition of bosons in one dimen- 
sion in such systems as Cooper pairs in Josephson junc- 
tions arrays^ and cold atoms in optical lattices^. Sine- 
Gordon model also provides a useful framework for un- 
derstanding properties of one dimensional spin systems 
with easy axis/plane anisotropics or in the presence of 
magnetic field (see Q for review). Interacting electron 
systems in one dimension often exhibit phase transitions 
between states with short range and power law corre- 
lations. Such transitions are also expected to be in the 
universality class of the quantum KT transition described 
by the quantum sine Gordon model 8 . Large classes of 
boundary-related phenomena 9 , including the problem of 
an impurity in an interacting electron gas (the so-called 
quantum impurity model) (see [l(| for review), Joseph- 
son junctions with dissipation^, Kondo-like models^, 
single electron transport in ultrasmall tunnel junctions 
have also been discussed in the framework of the bound- 
ary sine-Gordon model. Quantum sine-Gordon model 
provides a universal description for such wide range of 
systems because it is the simplest model with a gapped 
spectrum and relativistic low energy dispersion. 

One approach to understanding the sine-Gordon model 
is based on the renormalization group analysis (see 
Refs. [|[ and [HJ for the review). This method has been 
successfully applied to discuss thermodynamic and trans- 
port properties of two dimensional superfiuids and super- 
conductors^. 

Another more detailed approach is based on the in- 
tegrability of the model and existence of the complete 



exact solutio n 13 i 14 . From this analysis we know that the 
excitation spectrum of the model consists of solitons, an- 
tisolitons and their bound states-breathers. The number 
and properties of different breathers is controlled by the 
interaction strength (related to the parameter 7 intro- 
duced in section IB below). 

While theoretical understanding of the quantum SG 
model is quite advanced, very few experimental studies 
have been done that addressed dynamics of this system. 
Several experimental studies on one-dimensional magnets 
interpreted resonances in neutron scattering and ESR ex- 
periments as breather type excitation o 15 ' 16 ! 17 . All these 
examples include spin-1/2 magnetic chains perturbed by 
an effective <?-tensor and the Dzyaloshinskii-Morya in- 
teraction or by a staggered field in copper benzoat o 18 i 19 
and dimethylsulfoxideSS, and/or by strong magnetic field, 
like in copper pyrimidine dinitrateA 7 -. Apart from these 
measurements we do not have experimental evidences to 
verify our understanding of the spectrum of the quan- 
tum SG model. Considering the importance of the sine- 
Gordon model for understanding one dimensional quan- 
tum systems it is of great interest to find new, more direct 
approaches for experimental investigation of this funda- 
mental model. 

In this paper we show that a pair of coupled one- 
dimensional condensates provides realization of the quan- 
tum Sine-Gordon model. Such systems were recently 
realized in cold atoms in microtraps, where the RF 
potential controls the tunneling between the conden- 
sate s 21 ' 22 ' 23 . The advantage of using cold atoms to re- 
alize the quantum SG model is that cold atom systems 
are highly tunable. For example in Refs. [24ll26| it was 
demonstrated that the interaction and thus the Luttinger 
parameter describing one-dimensional condensates can 
be varied in a very wide range. Similarly the tunneling 
between the two condensates, which controls the gap in 
the sine-Gordon theory (see details below), can be tuned 
to a high precision. Coupled condensates can also be re- 
alized in optical lattices using superlattice potential s 21 ' 27 . 

One possibility to probe the excitation spectrum of the 
quantum sine-Gordon model is to study the response of 
the system to small periodic modulations of various pa- 
rameters of the model. In this paper we discuss three 



2 



types of such modulation experiments: modulation of 
the tunneling amplitude, modulation of the potential dif- 
ference between the condensates, and modulation of the 
phase of the tunneling amplitude (see Fig.[T]). By observ- 
ing the frequency dependence of the absorbed energy in 
various modulation experiments, one can measure struc- 
ture factors associated with appropriate perturbations. 
Here we theoretically compute structure factors for all 
three types of modulations (see Figs. [2l [4] [6] and [7]). As 
we demonstrate below, different types of modulations ex- 
pose different parts of the quantum sine-Gordon spec- 
trum and provide complimentary information about the 
system. We note that the calculations of structure fac- 
tors for several types of operators in the quantum sine- 
Gordon model was done in earlier papers in the context 
of condensed matter system a 28 i 29 ' 30 , mainly to describe 
properties of Mott insulating states. These results were 
summarized in a recent review 3 -^. Somewhat similar anal- 
ysis has been done in the context of the boundary sine- 
Gordon model (see Ref. [l(| for the review). In this paper 
we apply and extend this earlier analysis. 

We note in passing that in a separate publication we 
discussed a different possibility to study the structure 
of the sine-Gordon model using quench experiments 3 - 2 -. 
There we showed that the power spectrum of the phase 
oscillations after a sudden split of a single condensate 
into two contains detailed the information about various 
excitations in the sine-Gordon model. The quench and 
modulation experiments are thus complimentary ways to 
extract this information. 

This article is organized as follows. First we give an 
effective description of the system in Sec. Ill Al Then we 
briefly summarize general facts about the quantum sine- 
Gordon model in Sec. Ill Bl In Sec. IIIII we outline the 
explicit construction of the sine-Gordon form-factors. In 
Sec. IIVI we evaluate and analyze the structure factors 
for different kind of perturbation. Finally in Sec. [V] we 
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FIG. 1: Three types of experimental setups proposed to study 
linear response in two coupled condensates. The figure shows 
a projection of the confinement potential. In setup a) the 
tunneling amplitude is modulated: A(x,t) = A + Ai(a;,t) 
(see Eqs. (J5J, (|10p ): setup b) corresponds to the space- time- 
dependent population imbalance accessed via the change V(t) 
of the confinement depth in one of the well; in experiment of 
type c) the space-time modulation comes from winding of the 
phase of single particle tunneling, A(x,t) — Aexp(ik(t)x). 
A is parameter defined in Eq. The s-direction coincides 
with the longitudinal direction of condensates and is out of 
plane of the figure. 



discuss experimental implications and extensions of our 
work and summarize the results. 



II. MICROSCOPIC MODEL. 

A. Effective model of two coupled condensates 

We consider a system of two coupled interacting one- 
dimensional condensates. If interactions between atoms 
are short ranged, the Hamiltonian providing a micro- 
scopic description of the system is 
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(*\(x)y 2 (x) + t>Ux)*i(x)) , (1) 



where rij(x) — ^(x^j (x). Here t± is the coupling, 
which characterizes the tunneling strength between the 
two systems, and g is the interaction strength related to 
the 3D scattering length a^r> and the transverse confine- 
ment length a± = y // h/mu!± viai 3 - 3 - 
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where C — 1.4603 .... In what follows we assume that we 
are far from the confinement-induced resonance and the 
second term in the brackets in Eq. ((2|) is strictly positive. 
It is convenient to introduce a dimensionless parameter 
7, which characterizes the strength of interaction: 
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(3) 



where po is the mean field boson density. 

In the absence of tunneling, the Hamiltonian ([1]) corre- 
sponds to the Lieb-Liniger model 3 -^ of bosons with point- 
like interactions. For this model the effective Luttingcr 
liquid description 3 -^ agrees well with the exact solution 
for the ground-state properties 3 -^ and the low-energy ex- 
citations (see Ref. [331 f° r the review). We thus expect 
that the "bosonization" procedure underlying the Lut- 
tinger liquid formalism also provides a good description 
of coupled condensates, at least when the value of the 
tunneling t± is not too large. The advantage of this low 
energy description is that one can make explicit analytic 
calculations for both static and dynamic properties of the 
system. Physically, the main assumption justifying use of 
this effective Hamiltonian is the absence of strong density 
fluctuations in the system. At sufficiently short distances 
(shorter than the healing length) this condition is vio- 
lated. However, as one coarse grains the system, density 
fluctuations become weaker and the Luttinger liquid de- 
scription becomes justified. Thus the effective formalism 
can correctly describe phenomena at wavelengths longer 
than the healing length of the condensate. 
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Within the Luttinger liquid description the Hamilto- 
nian ([1]) splits into parts corresponding to individual con- 
densates Hi 2 and the tunneling between them H tU n- The 
first two are characterized by so called Luttinger liquid 
parameter K and the sound velocity v s , which we assume 
to be identical for both condensates: 



Hi: 



Hi 

2 



dx 



£ J I? i2 (a:) + — (0x01,2) 

A 7T 



(4) 



Here the phase field (f>x,2(x) is conjugate to the momen- 
tum Yli^(x). The relation between the original bosonic 

operators and the new fields is given by the bosoniza- 
tion rulei&2£ 



(5) 



where d<dj(x)/ir = p + Hj(x). 

In general, the relations of K and v s to the original mi- 
croscopic parameters m and g are not known analytically 
but can be deduced numerically from the exact solution 
of the Lieb-Liniger model. However, in the limit of large 
and small interactions one can obtain approximate ex- 
pressions^. Thus for 7 <C 1 
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i-^r l-r ( 6 ) 
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and in the opposite limit 7 1 



V F 1 - - 

7 



K 



1 + 1. 
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(7) 



Here «f = hnpo/m is the Fermi velocity. Note that 
v s K — vf as a consequence of the Galilean invariance. 
We also mention that the healing length £/, setting the 
scale of applicability of Eq. is approximately equal 
to the inter-particle distance: ^ ~ l//>o at 7 > 1 and 
becomes larger: £/, ~ 1/ajVt" at 7 <C 1. We note that ex- 
perimentally the coupling 7 can be tuned through a very 
wide ranged by either changing the density or transverse 
confinement length or by tuning a^o near the Feshbach 
resonance (see Eq. @). 

The tunneling part of the Hamiltonian H tun can be 
described effectively as 



-2A< 



2). 



(8) 



where we introduced the new parameter A. For small 
interactions 7 < 1 we have A rs txPo- At strong interac- 
tions A will be renormalized and can significantly deviate 
from its bare form. 

It is convenient to introduce the total and the relative 
phases <f>± = (</>i ± 2 )/v / 2 and n± = (Hi ± n 2 )/V2. 
Then the Hamiltonian splits into symmetric and anti- 



symmetric parts: H = H + + where 



H+ = y / <!.,■ 



^(n + (x)f + -(d x cj> + ) 2 

^(n_(,)) 2 + -(W-) 2 



7T 



2 A J dx cos(V20_) 



(9) 



Note that the symmetric part of the Hamiltonian H + in 
our approximation is completely decoupled from the anti- 
symmetric part H-. However, higher nonlinear terms in 
n and d x 4> will couple </>+ and 4>- so that the symmetric 
degrees of freedom can damp excitations of 0_ . 

The Hamiltonian H— can be simplified further. In par- 
ticular, we can rescale units of time: t — > tv s to set 
v s = 1. It is also convenient to rescale th e pha se field 
4> and the conjugate momentum: — ► Wtt/K and 
n_ -> ^if/Trn^. Then 



//. = ~ / das[II?.(s) + (0 K 0_) 2 - 4Acos(/30_)]. (10) 



where /3 = y/2-K/K. It is also convenient to introduce a 
parameter £: 



1 



(11) 



8tt - P 2 AK - 1 ' 
which we will use later in the text. 

B. Quantum SG model: general facts 

The spectrum of the quantum sine-Gordon Hamilto- 
nian (jTUJ) depends on the value of (3. For 47r < 1 < 8ir it 
consists of solitons and antisolitons. The point 1 = 4ir 
corresponds to the free massive fermion theory. In this 
case solitons and antisolitons correspond to particles and 
holes. For < (3 2 < 4ir the spectrum in addition 
to solitons and antisolitons contains their bound states 
called breathers. Note that within the Lieb-Liniger model 
K > 1, which implies that 1 < 2ir, so that we are al- 
ways in the latter regime^. The number of breathers 
depends on the interaction parameter £ and is equal to 
the integer part of l/£. We denote breathers by B n for 
n = 1 . . . [l/£]. In the gaussian limit of the sine-Gordon 
Hamiltonian (JTUJ) , where one expands cos(/30_) to the 
quadratic order in there is only one massive exci- 
tation corresponding to the lowest breather B\. Soli- 
tons (kinks), antisolitons (antikinks) and breathers are 
massive particle-like excitations. Soliton and antisoli- 
ton masses in terms of the parameters of the Hamilto- 
nian ((Tnj) were computed in Ref. [39(: 
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Breather masses are related to the soliton masses via 

M Bn = 2M S sin (j^j . (13) 

Note that at weak interactions (large K and small £) 
the lowest breather masses are approximately equidis- 
tant Mb„ oc n, which suggests that these masses corre- 
spond to eigen-energies of a harmonic theory. There is 
indeed a direct analogy between the breathers in the sine- 
Gordon model and the energy levels of a simple Joseph- 
son junction^. In the Josephson junction the energy 
levels also become approximately equidistant if the in- 
teraction (charging) energy is small. 

And finally at /3 2 = 87r the systems undergoes the 
Kosterlitz-Thouless transition so that for 1 > 8w the 
cosine term becomes irrelevant and system is described 
by the usual Luttinger liquid. 

C. Possible experimental probes 

We consider three different types of modulation exper- 
iments, which can reveal the spectrum of the quantum 
sine-Grodon model. We restrict ourselves to the linear 
response regime, which corresponds to small perturba- 
tions. Modulation of different parameters in the model 
focuses on different parts of qSG spectrum and provides 
complementary information. A schematic view of possi- 
ble experimental setups is given in Fig. [T] The direction 
of axes of two parallel condensates is out of the plane of 
figure. 

In this paper we discuss three types of periodic modu- 
lations: 

a) . Modulation of the magnitude of the tunneling am- 
plitude A(x,t). This kind of modulation couples to the 
potential density operator cos(/30). 

b) . Modulation of the relative potential difference be- 
tween the two wells V(x, t) . This modulation couples to 
the momentum operator II (x), which in turn is propor- 
tional to dt<j>. 

c) . Modulation of the phase of the tunneling ampli- 
tude, e.g. A*(x,t) — Aexp(±ik(t)x). After redefinition 
of variables <f> — > <f) — 4>±k(t)x/ (3 the term corresponding 
to the density of the topological current operator d x 4>- 
appears in the Hamiltonian. Note that the topological 
charge 

Q = iJ°° (14) 

is conserved and quantized as ±n with n being an integer 
(corresponding to presence of solitons and antisolitons). 
We note that in linear response the modulation of the 
type c) is equivalent to modulating the relative current 
between the two condensates. 

For each perturbation of the qSG model that we dis- 
cussed, one can define an appropriate correlation func- 
tion and susceptibility. According to the fluctuation- 
dissipation theorem energy absorption (per unit time) 



during modulation experiments is given by the imaginary 
part of these susceptibilities. The latter are also known as 
structure factors. In Sec. HVI we present typical structure 
factors corresponding to experiments of type a) - Fig. [21 
type b) - Fig.[H and type c) - Fig. [51 A general structure 
of the response function is given by a collection of peaks. 
These peaks can be classified into two groups: coherent, 
(5-function-like peaks coming from single breathers and 
the broader peaks coming from creating many (typically 
two)-particle excitations. From the brief overview given 
in Sec. IIIBI it is clear that the absorption spectrum sig- 
nificantly depends on the interaction parameter K . We 
give a detailed analysis of structure factors in Sec. HVl 

III. QUANTUM SINE-GORDON MODEL AND 
FORM-FACTORS OF ITS OPERATORS 

The Hilbert space of the sine-Gordon model can be 
constructed from the asymptotic scattering states. The 
latter can be obtained by the action of operators A ak (9) 
corresponding to elementary excitations on the vacuum 
state 

|0i0 2 ...0 n > ai)O2 ,... >On = <(0i)4 2 (^)...4„(e„)|o),(i5) 

where the operators A ak (6k) have internal index cor- 
responding to solitons (a*; = +), antisolitons {au = — ), 
or breathers (a^ = n, n = 1 . . . [l/£]) and depend on the 
rapidity 9^ which parametrizes the momentum and the 
energy of a single-quasi-particle excitation of mass M ak 
(see Eqs. (p j) . (fl3 |) ): E a = M a cosh(0), P a = M a sinh(0). 
These states are complete and therefore satisfy the com- 
pleteness relation: 




...0n|.(16) 



Elementary states or excitations in an integrable field 
theory can be described in terms of operators creating 
or annihilating asymptotic states. Their commutation 
relations involve the scattering matrix S 

A a (9 1 )A b (9 2 ) = S a J'{9 1 ~9 2 )A b ,{9 2 )A a ,{9 l ) 

At(0!)A+(0 2 ) = si b '{9 1 -9 2 )Al,{9 2 )Al{9 1 ) 

A a (9 1 )Al(9 2 ) = 2tt<WK0i - 2 ) 

+ S b b ' a U9 1 -9 2 )Al(9 2 )A al (9 1 ) (17) 

These relations are called the Zamolodchikov-Faddeev al- 
gebra. They generalize canonical commutation relations 
for bosons and fermions and reduce to them in some spe- 
cial cases. Since breathers are the bound states of solitons 
and antisolitons, the soliton-antisoliton (a = " + ",£> = 
" — ") scattering matrix has corresponding imaginary 
poles at 9 n = wr(l — n£), n = 1, ... [l/£]. 

Different components of S-matrix are related to sev- 
eral allowed scattering processes (related to the config- 
urations in the 6-vertex model) and as a building block 
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include the following quantity 

dx sm 



So(0) 



exp 



in (ff ) sinh (i^x 



sinh(x) cosh ( | 



(18) 



Form- factors 



The form factors F° of a given operator O of an 
integrable model are the matrix elements in asymp- 
totic states (p~5)) created by elements A and, of the 
Zamolodchikov-Faddeev algebra, Eqs. (fT7|) . Explicitly, 



F°(9 n . . . 0l)a„... ai = (0|O(0, 0)\9 n , . . . 6 1 ) a „ 



(19) 



Form-factors (FF) satisfy a set of axioms and functional 
equations, which together allow in principle to determine 
their explicit forms. Using the so-called crossing relation 
and the translation invariance all form-factors can be ex- 
pressed in terms of matrix elements of the form (I19[) : 

n 

{0\O{x,t)\e n ,...6 1 ) an ... ai =exp[iJ2iE j t + P j x)} 

x<O|0(O,O)|0 n ,...0i} OB ...« I , (20) 

where Ej = M 0j -cosh(%), Pj = M aj smh(9j) . Note 
that in the noninteracting limit of some generic model 
Zamolodchikov-Faddev algebra reduces either to Bose 
or Fermi canonical commutation relations. In this case 
form-factors simply reduce to the coefficients of the ex- 
pansion of the operator O in the second quantized form. 

Explicit analytical expressions for the form-factors de- 
pend on the specific type of the operator O. In this paper 
we will be interested in two particular types of operators: 
(i) the current operator in the sine-Gordon model: 



J" 



'2T 



(21) 



where /i = 0, 1 correspond to time t and space x compo- 
nents; and (ii) the other operator corresponding to the 
trace of the stress-energy tensor T^ v — d^d u (f> — S^C, 
where C is the Lagrangian, 



T = Tr(T{f) = 4Acos(/3<?!>_). 



(22) 



These two operators have different properties with re- 
spect to the Lorentz group and different charge parity. 
Therefore their properties can be distinguished by the 
topological charge (fl~4"]) and the charge conjugation oper- 
ator C defined by 

C|0) = |0), CAl^C- 1 = 4(d), (23) 

cAljetf- 1 = (-i)"4„(0) (24) 

The action of both the current operator and cos(/30_) do 
not change the value of Q. 



There is extensive literature on computation of form- 
factors for the quantum sine-Gordon model. To our 
knowledge such analysis was initiated in Ref. (40| and 
extended later in Ref . 1411 . These developments are sum- 
marized in the book [42l| and more recently in Refs. [43|, 

If there are bound states in the theory, like breathers 
in our case, the form-factors which include these bound 
states can be constructed using the residue at the poles of 
the soliton-antisoliton form-factors. To get form factors 
corresponding to higher breathers one can use a so-called 
fusion procedure. The receipt is that the breather B n+m 
appears as a bound state of the B n and B m breathers, 



B 



B,, 



B 



n+ra • 



Thus the residue of the pole in 
the form- factor FB n B m will correspond to the form- factor 
FB n+m ■ For example, to construct the form factor which 
includes the second breather B 2 we can use the form 
factor with n + 2 breathers B\ and the fusion formula 



iReSe^o^ai—anll ■ ■ • j Qm®n+l + iU\ 2 ~ — , 

6-n+i — iU\ 2 + -J = Ti 1 F ai .„a n 2(di, . . . , n ,6 n+ i), (25) 

where the fusion angle U\ 2 — 7r£/2 and the 3-particle 
coupling Y\ x = ^/2tan(7r£) is given by the residue of the 
S-matrix ([15]). 

1. Soliton-antisoliton form-factors 

1). O = e 1 ^ . For the sine-Gordon model different 
form factors for the phase operators e l °^ have been found 
in Refs. 42,43]. In particular, the nonvanishing "H — " 
(soliton-antisoliton) form factor is given by 



o W(0,0)i a t 



|4=(03)4:(0i)> = 0/jexp[J(0i 3 )] (26) 



2i cot 



(?) 



sinh ( 



112 



£ sinh 



'12+"T 



• exp =F 



?12 + 27T 



where 0\2 

1(8): 



9-i — On, and 



dt sinh 2 [i(l 



|sinh[t(£-l)] 



sinh[2t] cosh[i] sinh[t£] 



(27) 



Equations (|26p as well as the normalization of gen- 
eral cos(/3</>_) -type form-factors includes the following 
vacuum- vacuum amplitude: 46 



Q = (0|e^| 



(l+£)tan(f) T 2 



27rr2 i±« 



.(28) 



Higher order soliton-antisoliton form-factors (e.g. with 
two solitons and two antisolitons) are also available in the 
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literature^ 2 -. For generic values of K they are given by 
multiple integrals, whereas for the half-integer values of 
AK (in our notations) the expressions considerably sim- 
plify In the next section we explain that the relative 
contribution to the operator's expectation value from the 
states which include many soliton-antisoliton pairs is very 
small and can be neglected. 

2). O = dn,<p The form- factors for these operators 

can be deduced from the known form-factors of the cur- 
rent operator J M . In particular, the form- factor for the 
operator dt4>~ is 



Fl t < t-(9 1 ,e 2 ) = -F dt *-(9 1 ,6 2 ) 

4ir 3 / 2 V2M s cosh (gg£) sinh gj^) 
cosh (^|+hl) 

The FF for the d x 4>- is given by 



(29) 



F 



d t <t>- 



-F 



dt<l>- . 



4n 3 / 2 V2M s sinh 



(30) 



01+02/ 

2 / 



sinh ( 



01-02 N 







cosh 



24 



-1(0, 



Multi soliton-antisoliton form factors of this operator can 
be found in Ref. [42j in terms of integrals. Their relative 
contribution is small as well 2 ^ and we will neglect them. 



2. Breather form-factors 

1). O = e 4 ^. The one-breather form factors of the 
type (0\e ll3 ^\B n (0)} = Fs n can be computed from the 
residue of the soliton-antisoliton form-factors at points 
9 n = — n£)), or, equivalently, from the procedure of 
fusion of several breathers. For example the breather B 2 
is the bound state of two breathers B\ . This procedure 
is known as a bootstrap approach. 



F, 



cxp(-i/30)_ 



g f3 V2 co t(j ) sin(7mQ ex V [I(~9 n )}e^ 
cot(^)n' s r i 1 cot 2 (^) 



B„ 



(31) 

We note that Eq.(|3"T|) reveals the parity property of 
breathers: odd breathers are antisymmetric with respect 
to the charge symmetry transformation, whereas even 
breathers are symmetric. 

The form-factors of several breathers can be derived as 
well. Since the breather B\ is like a fundamental particle 
in a theory, all higher-level breathers form-factors can 
be expressed via a fusion procedure of Bi-form-factors. 
Because of the parity the non-zero form-factors for the 
operator cos(/3</>) contain either even breathers or even 
number of odd breathers. It is thus important to have 
an explicit expression for the B\ . . . i?i-form factors: 



Si. ..Si - 

- ^ 



\cos(p<p)\B 1 (6 2n )...B 1 (d 1 )) 
, det(S; j .) 
det (£„■)' 



(32) 



l<i<j<2n 



Here the (2n — 1) x (2n — 1) matrices 

sin[7r£(i — j + 1) 



S»j - <T2i-j, ■ - <T 2i -j- 



are expressed in terms of symmetric polynomials 

2n-l 



(2n-l) 



E 



il <i2< - <ik 



with Xi 



2 cos 



11(9) = Nexp 



2 



'2 sin 



2 



exp 



^ dt t 
1 1 2tt sin t 



dt sinh(t) sinh(t£) sinh(t(l + £)) 



sinh 2 (2i) 



x sinh I 1 — i— 

7T 



(33) 



N 



exp 



dt sinh(i) sinh(^) sinh(<(l + £)) 



sinh 2 (2t) 



Here the function 71(6) satisfies a useful relation which 
allows to evaluate its residue at 9 = — wr(l — n£) 



11(9)11(9 ± nr) 



sinh(6>) 



sinh(#) =F i sinh(7r£) 



(34) 



The formula (|32p can be shown to agree both with 
the sine-Gordon "bosonization" by Lukyano\*^i and the 
Bethe-Ansatz-based method of Ref. 45] . 

The form-factor for the B 2 — B 2 breathers can be de- 
rived from Eq. (|32p using the fusion procedure: 



F\ 



B2B2 



(9) 



1 



tan(7r£) 1 H , 

2 \ cosh(0) + cos«) 

1Z 2 (9)1Z(-9 - in£)K(-6 + iirQ 



2). O = dt4>-- Because of the charge reflec- 
tion symmetry, odd-type one-breather form-factors ex- 
ist for this operator. They can be found either di- 
rectly from the poles of the soliton-antisoliton form- 
factor ([H]) or via the relation (Q\d t <j)-\B n (9)) = 
lim a ^ ~M Bn cosh(6»)(0|e ia<?i - \B n (9)). The result is 



FTAO) 



n 



cot 

\ ■> I 

s=l 

M Bm V2(2^) 3 / 2 g 





tan 



7TS£; 



1/2 



(35) 



exp(J(-0 m )) cosh(0), 
1/2 



cot 

\ ■> i 

s=l 

M 5m V2(2^) 3 / 2 e 





(36) 



exp(/(-0 m )) sinh(0), 
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where m is odd. 

To compute form-factors for several breathers we note 
that there is a correspondence between the lowest B\ 
breather sector of the sine- Gordon theory and the sink- 
Gordon theory. Using this correspondence the B\ - form- 
factors can be obtained via the recursion formula^ from 
the form-factors of (^ s inh - operator in the sm/i-Gordon 
model. This recursion gives us odd-number non-zero B\ 
form-factors 



F 



dt4>- 



(2^) 3 / 2 A 2 " +1 g 
/?sin(7r£) 



M. 



Bi 



'2n+l 



(2n+l) „ , 
Xcr 2n+1 P2n+l(Xl, 



i<j 



cosh(#„ 

n=l 

i) 



(37) 



where Xi = e ; and first polynomials P are given by^ 

Pz{xi,X2,X3) = 1 (38) 

P 5 (xi,...,x 5 ) = cr 2 cr 3 -4 cos 2 «)cr 5 (39) 

These formulas allow us to compute, for example, the 
B\ — £?2 form-factor, which is the first nontrivial 2- 
particle form-factor after the soliton-antisoliton form- 
factor: 



F 



Bi b 2 VM = (M Bl sinh(^) + M B2 sinh(0 2 )) 



x F. 



F^X(Si2 



B1B2 

^ 3 / 2 A 3 £ 



12 J 



cos(^) (cosh(#) + cos(^)) K (-«tt(1 + 0) 



It is remarkable that the further fusion of £?i and B% 
particles produces the B$ form-factor, which is identical 
to the result from Eq. I|35p. This agreement is a direct 
manifestation of the bootstrap^. 



IV. LINEAR RESPONSE THEORY FOR 
COUPLED CONDENSATES 

Experimentally one can modulate the relative coupling 
between two condensates by applying small modulations 
in space and/or time. Moreover one can change the densi- 
ties in both condensates such that the relative density can 
vary. Following the previous discussion we will consider 
three relevant operators which describe different physical 
perturbations 



a) . Ot(x, t) = cos(/30_), 

b) . Oj 1 {x,t) = d t 4>-, 

c) . Oj (x,t) = d x (j)- 



(40) 
(41) 
(42) 



These operators directly correspond to the experimental 
setups (a), (b) and (c) discussed in the Introduction. The 



conserved topological charge in these notations is given 
by the integral of J°. 

In the linear response regime we are interested in com- 
puting two-point correlation functions 



(0(x,t)0t(O,O)), 



(43) 



where the expectation value is taken over the ground 
state of the unperturbed system. Inserting the complete- 
ness relation (|16p and using the relativistic invariance al- 
lows us to express this average in the following weighted 
sum over the intermediate states: 

r n Aei 

{0(I ,^(o,o)>=EE/n(^ 

xe ^U P^-B if |<o|O(o,O)|0„ . . . e 1 ) an ... ai \ 2 .(44) 

It is convenient to define the dynamical structure factor 
as a Fourier transform of this quantity, 

S°(q,uj) = Im[ / dxdte luJt - lqx (0{x,t)0^ (0,0))] 

J J — oo 

d6i 



(27r) n n! 



\F 



i / d\ ...a n | 



X 5 (q - M a 3 sinh 9 s jS\u-^2 M aj cosh 9j j , (45) 

where F° is the corresponding form-factor of the opera- 
tor O. We comment here on the structure of the expres- 
sion (|4"5)) . First, it is represented as a sum of contribu- 
tions coming from different excited states of the Hamilto- 
nian. Second, it is known from rather general phase-space 



argument: 



,48.4!) 



that the sum over intermediate states of 



the form (|45|) converges rapidly. In particular, it was 
shown in Ref. [29( that the two-soliton-two-antisoliton 
contribution is somewhat 500 smaller than the soliton- 
antisoliton one. Therefore, in general, the contributions 
with small number of excited particles will dominate the 
result . The expression for the structure factor ([43)1 can 
be also rewritten as the following sum 



S°(q,cj) 



Sf 1} {q,uj) + ^-Sf 2) {q,w) 



+ ... (46) 

The specific form of S^(q,u>) significantly depends on 
the charge parity of the operator and the topological 
charge parity. The current operator is odd with respect 
to the C-parity whereas cos(/3</>_) is even. All operators 
we consider do not change the topological charge. 

Explicit expressions for the S cos ^^ (q,uj) have already 
appeared in literatur o 30 i 31 in the context of ID spin-1/2 
systems discussed in the Introduction. We present them 
here for completeness. 

a). In this case O = cos(/?</>-). Then the single particle 
contribution to the structure factor is given by 



[l/£] 



cos(/30) 

s 2 „ 



S(s 2 



M B 2 



), (47) 
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S(q,o>) 


1.4 




1.2 




1 




0.8 




0.6 




0.4 




0.2 





masses M Al and Ma 2 can be generally expressed as 



FIG. 2: Structure factor for the experiment of type (a) cor- 
responding to the phase operator cos for the Luttinger 
parameter K — 1.15 and the tunneling gap A = 0.1. Peaks 
correspond to the following leading contributions (from left 
to right): a single breather B2; breathers B1-B1; and soliton- 
antisoliton A+A-. 



7 cos (j80) 




K 



of single-breather con- 



FIG. 3: Spectral weights \F C E 
tributions corresponding to the scheme a). Shown here are 
contributions for (from top to bottom) B2, -B4 and B§ (very 
small) . 



£(2) 



Re 



M 2 



Ml 2 



±MlM% 



where 



?12 = 9(q, ^) = arccosh 



- Ml - Mj 
2M Al M A . 2 



(49) 



(50) 



Direct substitution of the single breather form-factor, Eq. 
(f3"Tj) . the soliton-antisoliton form- factor (Eq.([2S])), and 
breather-breather form- factors (Eqs. (|3"2l . (135)) ) from the 
previous section gives the results for the corresponding 
structure factors. Note that the soliton-antisoliton con- 
tribution has to be weighted by a factor of 2. The re- 
sults are shown in Figs. [2] and |4] for some fixed values of 
the Luttinger parameter K and the inter-chain coupling 
A. We see that the generic form of the structure factor 
(and this of the absorbtion peaks) has a coherent part 
represented by a sequence of 5-peaks corresponding to 
coherent isolated even breather contributions and the in- 
coherent part corresponding to different two-particle con- 
tributions. In Figs. [2] and [4] we show the contributions 
from B\B\, A + A_, B2B2 (incoherent part) and B2, B4 
(coherent peaks) form-factors. 

S(q,w) 
1 



0.8 
0.6 
0.4 
0.2 



where s 2 — lu 2 — q 2 . Thus iSm corresponds to excitation 
of isolated even breathers. Here the spectral weights are 

^.cos(/30) _ |^,cos(/30) |2 ^g-j 



FIG. 4: Structure factor for the experiment of type a) cor- 
responding to the phase operator cos(/30) for K — 1.4 and 
tunneling A = 0.1. Peaks correspond to the following lead- 
ing contributions (from left to right): single breather B2; 
breathers Si Si; single breather B4,; and soliton-antisoliton 
A+A-. 



where F^^ ' is given in Eq. (f3"Tj) . So this contribution 
corresponds to the absorbtion peaks at id = oJB 2n (f° r 1 = 
0). The spectral weights are illustrated in Fig. ([3]) for n — 
1,2,3. Note that these weights decrease with increasing 
n as well as with increasing K. This plot suggests that 
breathers with small n dominate the absorbtion. 

The two particle contribution to the structure factor 
corresponding to excitation of particles Ai and A2 with 



b). For O = dt(f>- the general form of the structure 
factor is different: 

S*f(q,w) = Z%*_y5 (s 2 - M 2 ^) , (51) 

where Z^_ ± = \F°£_j l E B2n _,\ 2 is given by Eq. (03). 

The spectral weights Z^^ _ are compared in Fig. [5] for 
n = 1,2,3. 
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1 

0.75 
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0.25 



S(^,0.76) 

8 
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2 



K 



CO 



FIG. 5: Spectral weights \F^ ] \/E Bn u . , ; 
breather contributions appearing in setups b) and c). Shown 
here are contributions for (from top to bottom) Bi, B$ and 
B 5 . 



l^ tW |/-PsJofsingle- 



S(w,0) 

5 \ 



4 
3 
2 
1 



CO 



FIG. 7: Structure factor S(w, q = 0.76) for the setup scheme 
of type (b) corresponding to operator O = dt4> at K = 1.6 
and A = 0.05. The structure of peaks is the same as in Fig. [6] 
except for the appearance of the satellite breathers peak at 
small to. It starts from nonzero q and exists further for all g's. 



In Fig. [5] we show the structure factor for the oper- 
ator O = dt<p- for K = 1.6. As before we include contri- 
butions corresponding to excitations of B±, B3, A + A- 1 
and B1B2. 

c). The structure factor for this type of modulation is 
very much the same as for the case b). In particular, the 
w-dependence in the nominator of Eqs. ()5H53|) should be 
replaced by q. 



SUMMARY AND OUTLOOK 



FIG. 6: Structure factor S(w,q — 0) for the setup scheme 
of type (b) corresponding to operator O = dt<j> at K = 1.6 
and A = 0.05. The 5-peaks (solid bold line) corresponding 
to the breathers Bi (thick line) and B3 (thin line) are coher- 
ent contributions whereas peaks in the incoherent background 
(from left to right) correspond to the breather Bi — B2 and 
soliton-antisoliton contributions. 



Similarly the soliton-antisoliton contributions are 

I^12)| 2 



^Vs 2 -4M S 2 



cosh ( H|a 



cos [j 
(52) 

where #12 is defined in Eq. (f50|) with Ma x = Ma 2 = 
M s and 1(9) is given by Eq. ([2T|). And finally the two 
breather contributions are 



not 



(2),S 1 S 2 (9' a; ) 



Re 



u 2 \F(0 12 )\ 2 



2 _ «f 2 



Ml 2 ? 



,(53) 



{2M Bl M B2 y_ 

where #12 is defined in (|50p and F(0 12) is introduced in 



We considered a system of two coupled one dimen- 
sional condensates. We showed that the low-energy dy- 
namics of this system can be described by the quantum 
sine-Gordon model (Sec. Ill A)) . This model is integrable 
and supports collective excitations (solitons, anti-solitons 
and breathers). To reveal these excitations we propose to 
study modulations of the tunneling amplitude - type a), 
of the population imbalance - type b) and of the tunnel- 
ing phase - type c). Corresponding experiments provide 
complementary information about the structure of exci- 
tations of the quantum sine-Gordon model. The modu- 
lations of type a) reveal even-number coherent breather 
peaks and even two-breather contributions whereas mod- 
ulations of type b) and type c) couple to the odd sec- 
tor of the spectrum, so that the coherent peaks corre- 
spond to odd breathers and odd two-particle excitations. 
All types of modulations show the soliton-antisoliton re- 
sponse. The effect of the Luttinger interaction parameter 
K of the individual condensates is two-fold: (i) the spec- 
trum content is entirely determined by the strength of K, 
as it is explained in Sec. Ill Bl and (ii) the spectral weights 
of coherent single-particle excitations B n with n > 2 de- 
crease significantly with increasing K . This result is not 
surprising since at large K we anticipate that the low- 
est energy excitations are well described by the Gaus- 
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sian model of the massive scalar field. In this limit there 
are only massive phonon-like excitations corresponding 
to Bi. For the type a) modulation these excitations can 
be created only in pairs giving raise to the continuum 
contribution to the spectral function. While for type 
b) and type c) modulations the lowest energy contribu- 
tion comes from excitations of isolated B\ breathers and 
then there is a three-particle continuum. We note that 
in the weakly interacting limit K ^> 1 the isolated B\ 
peak corresponds to exciting Josephson oscillations be- 
tween the two condensates. Also we note that the soliton- 
antisoliton contribution to the spectral function rapidly 
diminishes with increasing K . . The experimental visi- 
bility of various contribution will also be determined by 
the tunneling strength A. In general one can observe an 
approximate scaling of structure factors with increasing 
A. 

We comment that the idea behind modulation exper- 
iments is similar to exciting parametric resonance in a 
usual harmonic oscillator. This analogy becomes even 
more transparent for the special case of zero dimensional 
condensates, which is equivalent to the Josepphson junc- 
tion. At weak nonlinearity the Josephson junction in turn 
is equivalent to a harmonic oscillator. The type a) mod- 
ulation of the tunneling amplitude is analogous to the 
modulation of the mass of this oscillator and the type b) 
modulation is similar to the modulation of the equilib- 
rium position of this oscillator. The strongest parametric 
excitation of a harmonic oscillator with frequency ujq oc- 
curs at u> = 2wo for the modulation of type a) and at 
u = uiq for the modulation of the type b). These transi- 
tions, in turn, correspond to excitations of the second -a) 
and the first - b) energy levels in the Josephson junction. 
As we discussed in Rcf. [32[ the first energy level of the 
Josephson junction corresponds to the B\ breather in the 
sine-Gordon model, the second energy level does to the 
B2 breather and so on. And indeed as we showed above 
the contributions from isolated B2 and B\ breathers to 
the absorbtion are dominant for the type a) and type b) 
modulations respectively. We can also come to similar 
conclusions using the classical analysi o 50 i 51 of the per- 
turbed SG model for small £ limit. One needs to bear 
in mind that here we deal with a nonlinear system and 
the parametric resonance strictly applies to a harmonic 
oscillator. So the response of the system to modulations 
is more complicated. However, qualitatively the picture 
remains very similar in this limit. 

Experimentally the absorbtion can be enhanced by in- 
creasing the magnitude of the modulation signal. Even 



though our analysis is limited only to weak perturbations, 
where one can use the linear response, we expect that the 
overall picture will not significantly change in the nonlin- 
ear regime. However, one always needs to make sure that 
the nonlinear effects do not cause various instabilities in 
the system. For example, the modulation of type c) is 
limited by the possible commensurate-incommensurate 
phase transition 5 -, which occurs if the breather's gap gets 
comparable with the change of the chemical potential 
per particle. Another effect which exists on the classical 
level is a dissociation of breathers into decoupled soliton- 
antisoliton pair for relatively strong perturbations 52 . 

It is important to realize that since we deal with an in- 
tegrable or nearly integrable system the energy absorb- 
tion is not necessarily related to the loss of the phase 
coherence, which is usually measured in standard time of 
flight experiments. Indeed, the two quantities are related 
if the absorbed energy is quickly redistributed among all 
possible degrees of freedom. However, in integrable or 
nearly integrable model this is not the case. For exam- 
ple, in a recent experiment by Kinoshita et. al^ it was 
demonstrated that there is no thermalization in a sin- 
gle one-dimensional Bose gas even after extremely long 
waiting times. Thus one possibility to measure the ab- 
sorbtion is to slowly drive the excited system into the 
non-integrable regime where the thermalization occurs 
and then perform conventional measurements. The other 
possibility is to directly measure the susceptibilities. For 
example, for the type a) modulation one can measure 
the expectation value of the phase difference between the 
condensates. For the type b) modulation one can mea- 
sure the population imbalance between the two conden- 
sates and for the type c) modulation one can measure the 
relative current between the two condensates. In princi- 
ple, one can always measure the loss of the phase contrast 
as it was done in Ref . [53| , however the sensitivity of such 
probe in integrable systems can be very small. 
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